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Abstract 

An accurate knowledge of the per-unit length impedance of power cables is necessary to correctly 
predict electromagnetic transients in power systems. In particular, skin, proximity, and ground return 
effects must be properly estimated. In many applications, the medium that surrounds the cable is not 
uniform and can consist of multiple layers of different conductivity, such as dry and wet soil, water, or 
air. We introduce a multilayer ground model for the recently-proposed MoM-SO method, suitable to 
accurately predict ground return effects in such scenarios. The proposed technique precisely accounts for 
skin, proximity, ground and tunnel effects, and is applicable to a variety of cable configurations, including 
underground and submarine cables. Numerical results show that the proposed method is more accurate 
than analytic formulas typically employed for transient analyses, and delivers an accuracy comparable 
to the finite element method (FEM). With respect to FEM, however, MoM-SO is over 1000 times faster, 
and can calculate the impedance of a submarine cable inside a three-layer medium in 0.10 s per frequency 
point. 


1 Introduction 

Electromagnetic transients, commonly induced by phenomena such as lightning and breakers operation, are 
a major source of power failures in today’s power systems [Tj . In order to understand and mitigate these 
transients, power engineers rely upon electromagnetic transient (EMT) simulation tools. EMT tools need 
broadband models for all network components including cables, which are increasingly used by the power 
industry. Broadband cable models MU can be obtained only if the per-unit-length (p.u.l.) impedance 
and admittance of the cable are accurately known. In particular, such parameters must precisely account 
for the frequency-dependent behaviour of the cable caused by skin, proximity and ground effects. 

Most EMT tools use analytic formulas M to compute cable impedance. Unfortunately, such formulas 
neglect proximity effects, which are strong in cables because of the tight spacing between conductors. Nu¬ 
merical tools based on the finite-element method (FEM) [7J [8] 0} [ED] or conductor partitioning [Ill |T21 [T5HTi] 
capture proximity effects inside the cable. However, these techniques require a fine discretization of the cable 
cross section in order to accurately capture skin effect, which leads to a large number of unknowns and long 
computational time. 

In various fault scenarios, current flows through the ground that surrounds the cable. Often, the ground 
surrounding the cable is non-uniform, and can only be modelled using multiple layers. For example, sub¬ 
marine cables are surrounded by layers of air, sea, and seabed, all with different material properties. In 
such scenario, the current flows inside each layer depending on each layer’s conductivity. Therefore, it is 
important to include effects of multilayer ground in cable models. 

*This work was supported in part by the KPN project ’’Electromagnetic transients in future power systems” (ref. 
207160/E20) financed by the Norwegian Research Council (RENERGI programme) and by a consortium of industry part¬ 
ners led by SINTEF Energy Research: DONG Energy, EdF, EirGrid, Hafslund Nett, National Grid, Nexans Norway, RTE, 
Siemens Wind Power, Statnett, Statkraft, and Vestas Wind Systems. 

tU. R. Patel and P. Triverio are with the Edward S. Rogers Sr. Department of Electrical and Computer Engineering, 
University of Toronto, Toronto, M5S 3G4 Canada (email: utkarsh.patel@mail.utoronto.ca, piero.triverio@utoronto.ca). 
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In most EMT tools, ground effects in underground cables are included through an approximation [T5] of 
Pollaczek’s formulas m- These formulas can only model a surrounding medium made up of two half-spaces, 
with the top layer being air, and the bottom layer being either soil or water. Pollaczek’s formulas cannot 
accurately capture effects of a multilayer ground. Several authors have developed analytic multilayer ground 
models mmm- These formulas must be used in conjunction with either analytic formulas or FEM to 
account for cables’ internal impedance j5]. FEM-based techniques |20| can model multilayer ground, while 
including skin and proximity effects. However, their computational cost is very high because a very large 
region around the cable must be meshed in order to accurately estimate ground return impedance. 

In pH [22], we presented an efficient technique called MoM-SO that computes proximity-aware p.u.l. 
impedance parameters of cables made up of round solid and hollow (tubular) conductors. Our technique was 
later refined in J23] to include the effects of single-layer ground and tunnels. The efficiency of MoM-SO stems 
from the use of a surface admittance representation of the cable. The surface admittance operator of the 
cable allows one to model the entire cable with a single equivalent current distribution on the boundary of 
the cable. The advantage of this is that we do not need to discretize the cross-section of the cable. Instead, 
only the edges of the cable are discretized, which leads to significant computational savings. With MoM-SO, 
ground effects are accounted for through the Green’s function of the surrounding medium. In this paper, we 
extend the technique to include a flexible multilayer ground model where the user can specify an arbitrary 
number of soil layers with different conductivity. Numerical results demonstrate the superior accuracy of 
the proposed model, and how it can simplify the task of modeling a power cable in a non-uniform soil. 

This paper is organized as follows. Section [2] defines the problem and the notation. In Sec. [3] we review 
the surface admittance formulation presented in [23] ■ Then, in Sec. [4] multilayer ground effects are included 
via a multilayer Green’s function. The final expressions for p.u.l. impedance inclusive of skin, proximity, and 
ground effects are provided in Sec. [5] Section [6] validates the proposed technique against FEM, and shows 
the importance of an accurate modeling of proximity and ground effects. 


2 Problem Definition 

Our goal is to compute the p.u.l. impedance of underground and submarine power cables. In order to handle 
a large number of scenarios, we will devise a method applicable to cables with the following characteristics: 

• any number of conductors, either solid or hollow, in arbitrary position and with arbitrary conductivity. 
Hollow conductors can be used to compactly represent an entire sheath or armour. Alternatively, 
individual strands can be also represented with solid conductors, owing to the excellent scalability of 
the proposed method; 

• a surrounding medium made by an arbitrary number of horizontal layers with arbitrary conductivity, 
as shown in Fig. |T| This versatile ground model can be used to describe ground effects in a large 
number of scenarios. For example, in a shallow submarine cable, the electromagnetic field produced by 
certain transients propagates partially in air, partially in water, and partially in ground. Such scenario 
can be easily described in the adopted ground model using the top-most layer to represent air, a middle 
layer for water, and the bottom-most layer for seabed; 

• the presence of a hole or tunnel around the cable can be modelled; 

• multiple cable systems, possibly buried in different holes, can be also modelled in order to investigate 
mutual coupling. 

For the sake of clarity, to describe the proposed method, we consider a cable made up of P solid round 
conductors, all with the same conductivity a, permittivity e, and permeability /i. Hollow conductors and 
conductors with different properties can also be handled, as shown in [23]. Conductor p is centered at 
{x p ,y p ), and has outer radius a p . All P conductors are placed inside a tunnel of radius a centered at ( x,y ), 
as illustrated in Fig [T] The hole consists of a lossless material with permittivity e, and permeability fi. The 
surrounding medium is made up of L layers, as shown in Fig. [T] Layer l has conductivity cq, permittivity 
Ei, and permeability /q. The top- and bottom-most layers, which are semi-infinite in the y-direction, are 
denoted as layer 1 and L , respectively. The layer in which the cable resides is denoted as layer s. 
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Figure 1: Cross-section of a simple cable made up of two solid conductors inside a hole. The medium 
surrounding the hole is modelled as L horizontal layers with different conductivity, permittivity, and perme¬ 
ability. 
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Figure 2: Left panel: simplified problem after the application of the equivalence theorem to the conductors. 
Right panel: problem after the application of the equivalence theorem to the hole. 


We are interested in computing the P x P p.u.l. resistance R(w) and inductance L(w) matrices that 
appear in the Telegrapher’s equation 

™=-[R{u>)+jwL{w)]I, ( 1 ) 

where vectors V = [Vi V 2 ... Vp\ T and I = [7i I 2 ... Ip] T contain the potential V p and the current 
I p in each conductor, respectively. 

3 Surface Formulation for the Cable 

In order to compute the cable impedance, we first apply the surface admittance formulation introduced 
in MM- This formulation reduces the complexity of the problem significantly, since a complex cable 
configuration can be described using a single equivalent current distribution. In this section, the surface 
formulation is briefly reviewed, and more details can be found in eh 122 mi- 
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3.1 Surface Admittance for Round Conductors 


We first expand the electric field on the boundary of each conductor using a truncated Fourier series. For 
the p -th conductor, this expansion is 


N p 

E Z (9 P ) = Y 4 p) ^' n V 

n=—N v 


( 2 ) 


are the Fourier coefficients of the electric field. The azimuthal coordinate 9 p is used to trace the 


where En' 1 

boundary of conductor p, as shown in Fig. 0 In@, N p denotes the number of harmonics used to represent the 
electric field distribution. For most practical cases, N p = 4 is sufficient 
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field on all the conductors are collected into a global vector E = 

Next, we replace all conductors with the surrounding hole medium, and introduce an equivalent current 
density on their boundary in order to keep the electric field outside the conductors unchanged. This trans¬ 
formation, illustrated in the left panel of Fig. [2] is enabled by the equivalence theorem [26]. As for the 
electric field in the equivalent current for each conductor is expanded in a truncated Fourier series 
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for all conductors are collected into a column vector 

r(i) hi) T (P) r(P)] T 

J -N i ’ ’ * J Ni ’ ’ * J -N P ’ * * J N P 


The current coefficients Jn' 1 

J = 

are related to E by a surface admittance operator Y s 

J = Y S E. 

This compact relation, which can be derived analytically for both solid [241 f2T] and hollow [! 
conductors, is sufficient to completely describe the conductors’ influence on the cable impedance. 


. As shown in [24J [21] , the equivalent current coefficients J 
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round 


3.2 Surface Admittance for a Cable-Hole System 

The geometry of the problem can be further simplified by applying the equivalence theorem a second time 
to the hole boundary. We first introduce the Fourier expansion of the magnetic vector potential on the hole 
boundary 

N 

A z (6) = Y An^ n§ ■ (5) 

n=—N 

Using the equivalence theorem [26], we replace the entire hole with the surrounding medium, i.e. having the 
same material parameters of layer s. An equivalent current density J s (9) is introduced on the hole boundary, 
and expressed in Fourier series 

N 
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~ 27TCL ^ J ' 


Jn9 
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where N controls the number of harmonics used to represent the cable-hole system. The hole equivalent 
current J s (6) is found with the following relation [23] 


J = Y S A + TJ 


( 7 ) 


where J = 


— N 


J 


N 


and A = 


-AT 




From (7), we can recognize two contributions to 


the equivalent hole current. The first term in ([7| is analogous to the surface admittance operator of round 
conductors Q, and models an empty hole. The second term in Q accounts for the presence of the cable 
conductors inside the hole. The transformation matrix T, maps the equivalent conductor currents © onto 
the boundary of the hole. Analytic expressions for Y s and T can be found in 
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Figure 3: Equivalent transmission line model used to represent the multilayer surrounding medium. Each 
layer is modelled as a segment of transmission line. 


4 Multilayer Ground Model 

Using the equivalence theorem, we have restored the homogeneity of the problem in each layer, as shown in 
the right panel of Fig. [2] In this equivalent configuration, the Green’s function can be conveniently used to 
relate currents and electromagnetic fields, and determine the cable impedance. We relate the vector potential 
and equivalent current on the hole boundary through the magnetic vector potential integral equation (26j 

A x (0) = -ft B J UO')G g (v{a,d),r{a,e'j)ade\ (8) 

where the integral kernel G g (.,.) is the Green’s function of the multilayer medium shown in the right panel 
of Fig. [ 2 ] Position vector r(a, 6) traces the contour of the hole. 


4.1 Multilayer Green’s function 

The Green’s function G g is related to the z-oriented magnetic vector potential A z (x,y) due to a z-oriented 
point source placed at (x r ,y')- The Green’s function of a multilayer medium can be found from the nonho- 
mogeneous Helmholtz equation [27l 28 : 

V 2 G s ( x, y) + k 2 s G g (x, y) = S(x - x',y - y') (9) 


in layer s, where the d-function is centered at (x',y r ) and k s = ujy s y/ui£ s — ja s is the wave number inside 
the layer. In all other layers, the Green’s function satisfies the homogeneous Helmholtz equation 


V 2 G g (x,y) + k?G g (x,y) = 0, 


( 10 ) 


where ki = ujyiy/u)£i — jcri is the wavenumber inside layer l. 
transform with respect to x to obtain 


To solve Q and (10), we apply the Fourier 


£^Gg(p x ,y) - (p 2 x - kl)G g (p x ,y) = e?^6(y - y’) (11) 

in layer s, and 

^G g (p x ,y) - (Pi - kf)G g (p x ,y) = 0 (12) 

in layer l ^ s, where 


Gg(Px 


/ OO 

Gg(x,y)e 3f3xX dx. 
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(13) 


It can be shown Bung that solving ( |11[ ) and ( |12[ ) is equivalent to solving the equivalent transmission 

line (TL) circuit shown in Fig. U| In this TL model, each layer of the background medium is modelled as a 

'—' —1 /2 

segment of TL, with length equal to the height of the layer, characteristic impedance Zi = ( Pi — kf) , 

and propagation constant 7 1 = ^PI — kf. The point source in ^ is modelled as a current source I s = e^ xX ' 
located at y = y'^ In the TL model of Fig. [ 3 ] the voltage along the line is equal to the spectral domain 
Green’s function G g (P x ,y). 
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Figure 4: Simplified transmission line model assuming that both source coordinates (x',y') and observation 
coordinates (x,y) are in layer s. 


To evaluate ( |8|), we only require voltage in layer s of the TL model in Fig. [3j Therefore, we can simplify 
the model in F igT[3] by replacing layers s +1,..., L by an equivalent impedance Z eq ^ s+ i and layers 1,..., s — 1 
by an equivalent impedance Z eq ^ s - 1 , as shown in Fig. [4j Both equivalent impedances are easily calculated 
using TL input impedance formulas found in most electromagnetic and TL theory textbooks [29] . The 
solution of the equivalent circuit in Fig. [4] is 1501 
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where T ^ and Tr are the reflection coefficients 
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Finally, we take the inverse Fourier transform of (14) to obtain the desired Green’s function 
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G g (x,y) = — / G g (J3 x ,y)e-^ x d^ , 


where the integral can be evaluated numerically. 
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4.2 Discretized Integral Equation 


Next, we substitute the Fourier series expansion of vector potential ([5]) and equivalent hole current ([6]), as 
well as the multilayer Green’s function © into Q. The resulting integral equation is discretized using the 
method of moments ED, obtaining 


A = ~n ( 


GgJ , 


where G g is a matrix representation of the discretized multilayer Green’s function 


into (18), we obtain the magnetic vector potential coefficients 
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. By substituting ([T]) 
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Figure 5: System of three single core cables used for validation in Sec. [ 6 ] 


Table 1: Single core cables of Sec. [ 6 ] geometrical and material parameters 


Core 

Outer diameter = 39 mnr, p = 3.365 • 10 -s Q ■ m 

Insulation 

t = 18.25 nrm, £ r = 2.85 

Sheath 

t = 0.22 nrm, p = 1.718 • 10 -s • m 

Jacket 

t = 4.53 mnr, e r = 2.51 


5 Computation of Cable Impedance 

To compute the p.u.l. impedance parameters we need the electric field on the conductors, which can be 
found by evaluating the electric field integral equation [26 


E z {Y p (9p)) 


-jujA z 


dV 

Hz' 


( 20 ) 


on the contour of each conductor. In (20), the position vector r P (9 P ) traces the contour of the p-th conductor. 
Using the same steps outlined in |23| . we finally arrive at the p.u.l. resistance and inductance of the cable 


R(w) =Ke|(u T (l- JW Y s ^)- 1 Y s u) 'j , 

L(w) =w- 1 Im|(u T (l-iwY s ’I')- 1 Y s u)” 1 | , 


where 


’S' = HDi 


s G g Y, 
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g T — /tGc 


AG c, 


and matrices H, D 1; 1, U, G 0 , and G c are defined in [25] , 


( 21 ) 

( 22 ) 


(23) 


6 Numerical Examples 

We now present two numerical examples to validate the proposed technique and to emphasize the need to 
include proximity-aware cable and multilayer ground models in EMT simulations. 

6.1 Example #1 - Three Tightly-Spaced Single-Core Cables 

6.1.1 Cable Geometry and Material Parameters 

The first example consists of three single-core cables (SC) buried under a shallow sea, as shown in Fig. [5] 
The geometrical and material properties of the cables are listed in Table [T] The surrounding medium is 
modelled using L = 3 layers. The top layer models air. The second layer models the sea, with height of 
10 nr, conductivity cr 2 = 5 S/m and electrical permittivity £2 = 81 £ 0 , which are typical for sea water 2)21 . 
The bottom layer represents a seabed with conductivity <73 = 0.05 S/m and permittivity £3 = 15 £0 [32] , 
The three SC cables are placed 1 nr below the sea-seabed interface, as shown in Fig. [5} The center-to-center 
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Figure 6: Cable system of Sec. 6.1 zero-sequence inductance (top panel), and resistance (bottom panel) 
obtained with MoM-SO (o), FEM (•), and analytic formulas (— ). Cable screens are left open. 


distance between adjacent cables is D = 85 mm. Since the SC cables are touching each other, significant 
proximity effects between them are expected. 

6.1.2 Impedance Validation 




Figure 7: As in Fig. [6j but when a positive-sequence excitation is applied to the core conductors. 

We first calculated the 6x6 impedance matrix of the cable using the proposed MoM-SO approach. In 
MoM-SO, the discretization parameters N p and N were set to 4. Next, we repeated the computation with 
a commercial FEM solver (COMSOL Multiphysics [33]), following the approach presented in [20]. In order 
to reproduce skin effect adequately, boundary layer elements were used to finely mesh the conductor edges. 
In FEM, we also used the infinite element domain to truncate the surrounding medium. A total of 246,818 
elements were required in the FEM simulation to mesh the cable and the surrounding medium. Finally, we 
calculated the impedance of the cable with the analytic formulas (cable constants [5j) that are implemented 
in most EMT tools. 

In order to facilitate the comparison of the results obtained with MoM-SO, FEM, and analytic formulas, 
we reduced the 6x6 impedance matrix to a 3 x 3 matrix by assuming that the screens of the cables are 
open at both end^j i.e. we assumed that there is zero net current inside the screens. Figures [6] and [7] 
show the cable inductance and resistance obtained with MoM-SO, FEM, and analytic formulas when zero- 
and positive-sequence excitations are applied to the core conductors. An excellent agreement between the 
results obtained with the proposed MoM-SO method and FEM can be observed. Analytic formulas instead 
return inaccurate parameters in both cases. In the case of positive sequence excitation (Fig. 7|), inaccuracy 
is mainly attributed to the neglection of proximity effects. In the zero sequence case (Fig. re) inaccuracy is 
due to proximity effects and the lack of an accurate multilayer ground model. 

1 This assumption is made only to simplify the comparison of the cable parameters, and will not be used in the subsequent 
transient simulations. 


8 









Table 2: Example of Sec. 6.1 CPU time required to compute the impedance at one frequency* 


Test Case 

MoM-SO 

FEM 

Speed-up 

Three-layers 

0.105 

183 

1743 X 


*Simulations were run on a system with 16 GB memory 
and 3.40 GHz processor. 
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Figure 8: Cross-bonded cable system setup. 


6.1.3 Timing Comparison 

The simulation times to calculate the impedance of the cable are summarized in Table [2] FEM took 183 s 
per frequency to calculate the impedance of this cable. On the otherhand, MoM-SO took just 0.10 s per 
frequency to calculate the impedance parameters with the same accuracy. 

6.1.4 Transient Simulation - Crossbonded Cable 

Next, we compare the transient waveform predicted with analytic formulas and MoM-SO. We consider the 
setup in Fig. [8] where a unit step excitation is applied to the core conductor of the left-most SC cable 
in Fig. [5] We created two universal line models j4] for the cable. The first model was derived from the 
impedance obtained with MoM-SO, while the second model was derived from the impedance calculated with 
analytic formulas. For both models, shunt admittance was calculated using the formulas from |34j . Figure [9] 
shows the transient voltages at node 1 and node 2 predicted with MoM-SO and analytic formulas. The 
voltage waveforms predicted with analytic formulas significantly deviate from the waveforms predicted with 
MoM-SO, as a result of the neglection of proximity and multilayer ground effects. These results show how 
the proposed method can lead to more accurate transient results with respect to existing EMT tools, which 
are mostly based on analytic formulas. 

6.2 Example # 2: Three Widely Separated Single-Core Cables 

6.2.1 Geometry 

The first example showed the influence of proximity and multilayer ground effects on cable impedance and 
transient results. In this example, we increase the separation between SC cables to 2 m, thereby minimizing 
proximity effects. We compare three different ground models: 

• three-layer air-sea-seabed model, where the surrounding medium is modelled as in Fig. [5} 

• two-layer air-sea model, where the presence of the seabed is neglected, and the sea layer extends to 
V = - 00 ; 

• two-layer sea-seabed model, where the presence of air is neglected, and the sea layer extends to y = 00 . 

The two-layer models represent what is currently possible with most EMT tools, which model ground as a 
single medium or as a two-layer medium m- 


9 










































Figure 9: Cable System of Sec. 6.1 node 1 and node 


formulas ( — ) for the setup shown in Fig. [8] 


2 voltages obtained with MoM-SO (—) and analytic 




Figure 10: Cable system of Sec. 6.2 zero-sequence p.u.l. inductance (top panel) and resistance (bottom 
panel) obtained with the three-layer air-sea-seabed model in MoM-SO (o), three-layer air-sea-seabed model 

in FEM (•), two-layer air-sea model (- - ), and two-layer sea-seabed model ( - ). The screens are continuously 

grounded. 


6.2.2 Grounded Screens 


We calculated the impedance of the cable using MoM-SO with the three ground models. For validation, the 
computation was also performed with FEM in the three-layer case. The 6x6 matrix was then reduced to a 
3x3 impedance matrix by assuming that the screens of the cable were continuously grounded. Figures 10 
and EH show the p.u.l. inductance and p.u.l. resistance values obtained using both MoM-SO and FEM 
(three-layer model only) when zero- and positive-sequence excitations are applied to the core conductors. 
Figures [TU] and EH show an excellent agreement between the impedance calculated with FEM and MoM-SO 
using the three-layer air-sea-seabed model, which validates the proposed approach. On the other hand, we 
note that the two-layer models underestimate the zero-sequence inductance at low frequencies. These results 
confirm the superior accuracy of the proposed method with respect to existing techniques puna. 


6.2.3 Open Screens 

We now assume that the cable screens are open. Figures [12] and [13] show the zero- and positive-sequence in¬ 
ductance and resistance obtained by exciting core conductors with zero- and positive-sequence currents. The 
parameters obtained with MoM-SO match very well the reference FEM results. The two-layer air-sea model 
produces inadequate results because at low frequencies the zero-sequence inductance is underestimated, and 
at high frequencies the positive-sequence inductance is underestimated. The two-layer sea-seabed model is 
more accurate than the air-sea model, but still underestimates the zero-sequence inductance at low frequency. 
It should be noted that the two layer sea-seabed model cannot be utilized in some EMT tools that require 
the conductivity of the top layer to be zero. 
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Figure 11: As in Fig. 10 


but when a positive-sequence excitation is applied to the core conductors. 
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Figure 12: Cable system of Sec. 6.2 zero-sequence inductance (top panel), and resistance (bottom panel) 


obtained with the three-layer air-sea-seabed model in MoM-SO (o), three-layer air-sea-seabed model in FEM 
(•), two-layer air-sea model (- - ), and two-layer sea-seabed model (-). Cable screens are left open. 


6.2.4 Timing Comparison 

In MoM-SO, discretization parameters N and N were set to 4. In FEM, the cable and surrounding medium 
were meshed with 255,380 elements. The timing results for both FEM and MoM-SO simulations are sum¬ 
marized in Table [3] MoM-SO takes only 0.1 s per frequency, and is 1661 times faster than FEM, which 
requires almost three minutes per frequency. 


Table 3: Example of Sec. 6.2 


CPU time required to compute the impedance at one frequency* 


Test Case 

MoM-SO 

FEM 

Speed-up 

Three-layers 

0.108 

179 

1661 X 


*Simulations were run on a system with 16 GB memory 
and 3.40 GHz processor. 


6.2.5 Transient Simulation - Crossbonded cable 

Finally, we compute the transient voltages excited by a unit step voltage applied to the three SC cables 
in cross-bonded configuration (see Fig. [8|. Figure 14 shows the transient voltages at nodes 1 and 2. The 
results show that there is up to 40% deviation between the results obtained with the air-sea model and 
the air-sea-seabed model. In this case, the two-layer sea-seabed model returns accurate results. However, 
if the water depth is reduced, this model becomes inaccurate. Since no two-layer model is accurate under 
all possible cases, their use requires the EMT engineer to understand which one is more appropriate for a 
certain cable. With MoM-SO, this dilemma is avoided, and accurate cable parameters are computed in less 
than a second, making cable modeling a straightforward and less error-prone task. 
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Figure 13: As in Fig. [l2j but when 
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Figure 14: Example of Sec 6.2 voltages predicted at nodes 1 (top panel) and 2 (bottom panel) of the 


configuration in Fig. [8j Plots compare the results obtained with three different ground models: three-layer 
air-sea-seabed(—), two-layer air-sea (- - ), and two-layer sea-seabed model (o). 


7 Conclusion 

A multilayer ground model was proposed for the MoM-SO approach for cable impedance calculation. With 
the proposed model, the non-uniformity of the medium which surrounds submarine and underground cables 
can be accurately taken into account. Numerical results show that the proposed method leads to better 
transient predictions than analytic formulas currently used in most electromagnetic transient simulators. 
While the level of achievable accuracy is comparable to a finite elements analysis, MoM-SO is more than 
1000 times faster than finite elements. Moreover, it is simpler to use, since it is fully automated and avoids 
meshing-related issues. We believe that these improvements make accurate cable modeling a simpler task 
for both transients experts as well as power engineers in general. 
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